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Abstract 

We study the diffusion of heavy quarks in the Quark Gluon Plasma using the Langevin equations 
of motion and estimate the contribution of the transport peak to the Euclidean current-current 
correlator. We show that the Euclidean correlator is remarkably insensitive to the heavy quark 
diffusion coefficient and give a simple physical interpretation of this result using the free streaming 
Boltzmann equation. However if the diffusion coefficient is smaller than ~ 1/(ttT), as favored by 
RHIC phenomenology the transport contribution should be visible in the Euclidean correlator. 
We outline a procedure to isolate this contribution. 
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I. INTRODUCTION 



The experimental relativistic heavy ion program has produced a variety of evidences which 
suggest that a Quark Gluon Plasma (QGP) has been formed at the Relativistic Heavy Ion 
Collider (RHIC) [1, 2]. One of the most exciting results from RHIC so far is the large az- 
imuthal anisotropy of light hadrons with respect to the reaction plane, known as elliptic flow. 
The observed elliptic flow is significantly larger than was expected from kinetic calculations 
[3], but in fairly good agreement with simulations based upon ideal hydrodynamics [4-8]. 
This result suggests that the transport mean free path is small enough to employ thermody- 
namics and hydrodynamics to describe the heavy ion reaction. However, this interpretation 
of the RHIC results demands further theoretical and experimental corroboration. 

Experimentally, this interpretation can be challenged by measuring the elliptic flow of 
charm and bottom mesons [9-11]. The first experimental results show a non-zero elliptic 
flow for these heavy mesons. Naively, since the quark mass is significantly larger than the 
temperature of the medium, the relaxation time of heavy mesons is ~ M/T longer than the 
light hadron relaxation time 

heavy light 

Consequently the heavy meson elliptic flow should be reduced relative to the light hadrons. 
Recently, a variety of phenomenological models have estimated how the transport mean 
free path of heavy quarks in the medium is ultimately reflected in the elliptic flow [12-14]. 
The result of these model studies is best expressed in terms of the heavy quark diffusion 
coefficient. (In a relaxation time approximation the diffusion coefficient is related to the 
equilibration time, r]j eavy = j^D.) There is a consensus from the models that if the diffusion 
coefficient of the heavy quark is greater than 

the heavy quark elliptic flow will be small and probably in contradiction with current data. 

Theoretically, transport coefficients have been computed in the perturbative quark gluon 
plasma using kinetic theory [15, 16]. The heavy quark diffusion coefficient has also been 
computed [12, 17, 18]. Recent efforts have also explored some meson resonance models and 
found a substantially smaller diffusion coefficient than in perturbation theory [19]. The 
ambiguity in these calculations underscores the need for reliable non-perturbative estimates 
of transport coefficients in the QGP. 

Kubo formulas relate hydrodynamic transport coefficients to the small frequency behavior 
of real time correlation functions [20, 21]. Correlation functions in real time are in turn 
related to correlation functions in imaginary time by analytic continuation. Karsch and 
Wyld [22] first attempted to use this connection to extract the shear viscosity of QCD 
from the lattice. More recently, additional attempts to extract the shear viscosity [23, 24] 
and electric conductivity [25] have been made. We will argue on general grounds that 
Euclidean correlations functions are remarkably insensitive to transport coefficients. For 
weakly coupled field theories this has been discussed by Aarts and Resco [26]. For this reason, 
only precise lattice data and a comprehensive understanding of the different contributions 
to the Euclidean correlator can constrain the transport coefficients. 
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In this paper we are going to estimate the contribution of heavy quark diffusion to Eu- 
clidean vector current correlators. The case of heavy quarks is special since the time scale 
for diffusion, M/T 2 , is much longer than any other time scale in the problem. In terms of the 
spectral functions, this separation means that transport processes contribute at small en- 
ergy, iq ~ T 2 /M, and all other contributions (e.g. resonances and continuum contributions) 
start at high energy, u > 2M. For light quarks, transport contributes to meson spectral 
functions for u ~ g 4 T. This scale is separated from the energy scale of other contributions, 
uj ~ T, gT, only in the weak coupling limit g <^1. 

The behavior of vector current correlators at large times can be related to the heavy 
quark diffusion constant. Euclidean heavy meson correlators at temperatures above the 
deconfinement temperature have been calculated on the lattice and attempts to extract 
spectral functions have been made [27-29]. Transport should show up as a peak at very 
small frequencies, u ~ 0. So far, it has not been observed in these studies. Obviously, it is 
very difficult to reconstruct the spectral functions from the finite temperature lattice corre- 
lators, as the time extent is limited by the inverse temperature. However, the temperature 
dependence of the correlators can be determined to very high accuracy [29, 30] and therefore 
some information about the transport can be ascertained. 



II. LINEAR RESPONSE AND THE SPECTRAL DENSITY 

This section briefly reviews linear response which is the appropriate framework to connect 
the Langevin and diffusion equations to the current-current correlator [20]. We will also 
define the spectral density which is needed to relate the Euclidean current-current correlator 
measured on the lattice to its Minkowski counterpart. 

Consider a small perturbing Hamiltonian 

H = H - J d 3 x h(x, t) 0(x, t) , (2.1) 

where /i(x, t) is a classical source. Now imagine that we slowly turn on the external source 
h(x,t), and then abruptly turn it off at time t — 0. h(x.,t) obeys 

h(jc,t) = e et 6{-t)h°{x). (2.2) 

The expectation value of (SO(x.,t)) in the presence of the perturbing Hamiltonian is 

(*0(x,f)> = +iJ d 3 y£ dt' ([O(x,i),O(y,0]> Ky,t'). (2.3) 
Using translational invariance and taking spatial Fourier transforms we have 

/+oo 
dt' X (Kt-t')h(k,t'), (2.4) 
-oo 

where 

x(M-f) = / d 3 xe-^z6(t-t') ([O(x,f),O(y,0]> , (2-5) 
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is the retarded correlator. When confusion can not arise we use momentum labels p, k, q, . . . 
rather than position labels x, y, z, . . . to distinguish the spatial Fourier transform of a field 
(O(M)) = JV kx (0(x,t)) from the field itself, (0(x,t)). 
For t > 0, differentiating with respect to £ we have 

.9 /" + °° d 

- (SO(k, £)> = J dt' - X (k, £ - £') h(k, £') . (2.6) 

Using J^x(k, £ — £') = — ^x(k, t — £'), integrating by parts with respect to £', and using 
Eq. (2.2), we find a relation between expectation values and correlators 

^(50(k,t)) = - x (k,t)h°(k). (2.7) 

The external field h°(k) can be eliminated by using the relation between the static suscep- 
tibility Xs, the initial condition (50(k,t)), and the external field 

(5O(k,t = 0)}=Xs(k)h°(k), (2.8) 

where the static susceptibility Xs(k), follows from Eq. (2.5) 

POO 

Xs (k)= / dt'e-^'xiKt'). (2.9) 
Jo 

Eliminating the field h°(k), we find 

X s (k)^- t (SO(k,t)) = - X (k,t) (50(k,t = 0)) . (2.10) 

This result relates the time evolution of an average from a specified initial condition to an 
equilibrium correlator %(k, i). 

The function x(k, t) is related to the spectral density. The Fourier transform of the 
retarded correlator can be written 

r+oo 

X(k,u)= / dte + ^ X (k,t). (2.11) 
Jo 

x(x, t) is real, and since the integration is only over positive times, ou) is analytic in 
the upper half plane. Provided the Hamiltonian is time-reversal invariant and the operator 
O has definite signature under time reversal, ([0(x, t) ,O(y,0)]) is an odd function of time 
and x(k, t) is an even (odd) function of k (time). The spectral density, p(k,u), is defined 
as the imaginary part by n of the retarded correlator 

p(k,u) = ImX ^' a;) =^Jd 3 * e-* k - x +^([0(x,t), 0(0,0)]) . (2.12) 

By inserting complete sets of states, one may show that the spectral density is an odd 
function of frequency and is positive for u > [31]. 

The Euclidean correlator may be deduced from the spectral density. Euclidean tensors are 
defined from their Minkowski counterparts, 0$"*" v v (-it) = (-i) r (i) s v v (r), 
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where r and s are the number of zeros in {/ii . . . //„} and {u\ . . . u n } respectively. In what 
follows, we will drop the "M" on Minkowski operators but indicate on Euclidean op- 
erators. With these definitions and Euclidean tensors transform under 
0(4) in the zero temperature limit. Correlators in Euclidean space time are of the following 
form: 

G(k,r) = J rf 3 xe ,k ' x (O e (x,t)O e (0,0)) = {-l) r+s J rf 3 xe lk ' x J D > (x, -it) , (2.13) 

where -D > (x, t) = (0(x, t)O(0, 0)) . Usually, the lattice works with at zero spatial momentum 
k = 0. In Minkowski space, we work with the Fourier transform of -D > (x, t), 

D > (k,cu) = J d i xe +wt -^D > {^t). (2.14) 

Similarly, we define -D < (x, t) = (0(0,0) 0(x, t)) and its Fourier transform. Thus, the spec- 
tral density, Eq. (2.12), is given by 

P (k,,) = (215) 

Using the Kubo-Martin Schwinger (KMS) relation -D > (k, t) = D < (k, t+i/T), and its Fourier 
counter-part -D > (k, oo) = e +u) l T D < (k ) uj), one discovers the relation between the spectral 
density and the Euclidean correlator, 

G (k, T) = ( - iy+ .jy p(k , w) ^M^M . (2 , 6) 

Again, given an operator, O^ 1 '"^ Vl ... Vn , r and s are the number of zeros in the space-time 
indices {/ii . . . fi n } and \v\ . . . u n } respectively. 

For our discussion two correlators will be important: the density-density correlator 

D> JV (x,t) = <J°(x,t)J°(0,0)> , (2.17) 
and the current-current correlator 

D>^(x,t) = <J 4 (x,t)^(0,0)> . (2.18) 
These correspond to the Euclidean correlators calculated on the lattice 

GW(x, r) = < J° E (x, r) J° (0, 0)> = -D> N (x, -it) , (2.19) 

G&(x, r) = (4(x, r)4(0, 0)> = £>>«(x, -ir) . (2.20) 

The corresponding retarded correlators Xiviv( x , and Xjj( x >0 can be introduced in the 
same way. The Fourier transforms of current-current correlators can be decomposed into 
longitudinal and transverse parts. For the retarded correlator we write: 

/ p Uj \ Li Uj 

xU^ = ( -p- ~ ^ ) *-( k > <") + "ifcT ^-( k ' <") • ( 2 - 21 ) 
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Current conservation relates the density- density and the longitudinal current-current corre- 
lators 

u 2 k % W 

— XNN (k, u) = -^rxKK u) = XjAK • (2.22) 

For k = there is no distinction between the longitudinal and transverse parts and therefore 
for k <C T , Xjj(k, w) ~ xjj(k, w). Since the transverse component of the current-current 
correlator is not studied in this work, we will drop the "L", and for instance, G 3J and p, 7J 
are short for Gj d and p L ir 

At finite temperature the spectral function can be written as 

A/J (k, u) = p 1 / 7(k, co) + pff(k, u) , (2.23) 

where the last term is just the zero temperature part and the first term is the low energy 
lu <C T 2 /M contribution. In the next two sections we will discuss how to estimate the low 
frequency part. 



III. TRANSPORT IN EUCLIDEAN CORRELATORS 

In this section we estimate how the low frequency part of the spectral function contributes 
to the Euclidean current-current correlator. To leading order, the moments of the spectral 
function, the time derivatives of the retarded correlator at t — 0, and the derivatives of the 
Euclidean correlator at (3/2, are in one-to-one correspondence. Since the leading contribution 
is due to time derivatives at t — 0, the free streaming of heavy quarks gives the dominant 
contribution to the Euclidean current-current correlator. The first scattering correction 
appears in the second (fourth) derivative at (3/2 of the current-current (density-density) 
euclidean correlator. 

First let us start with the density- density correlator. For small frequencies, the kernel in 
Eq. (2.16) is given by 2T/u, and thus we can write 

-G l Z(k,r)^2T r — p l Z(k,u;). (3.1) 
Jo u 

Inserting the definition of the spectral density 

1 f°° 

ff%(k,u) = - dt sin(^) W (k,t) , (3.2) 

71" Jo 

and performing the integral over frequency, we find 

POO 

-G l Z(k,r) ~ / dt X l Z(Kt) = T Xs (k) . (3.3) 
Jo 

The last equality follows from Eq. (2.9). Similarly, the low energy contribution to the 
longitudinal current correlator is 



f°° duj 

G 1 7 °7(k,r)^2T / — p 1 7 T(k,^) 
Jo u 



(3.4) 



6 



Inserting the spectral density 



k,u;) = 



nk- jo 

and performing the integral over frequency, we find 



dt sm(wt) x NN (Kt) , 



(3.5) 



G57(k,r) = ^ 



^ (1) X ^(k,t) + 



1 



24 T 2 



a t at 



,(k,t) 



cf) w (k,t)-(r- / 3/2) 2 + ... 



t=o 
(3.6) 

Thus we see that the dominant low energy contributions to the Euclidean correlator is given 
by the short time behavior of the retarded correlator. Indeed, as seen from Eq. (3.4) and 
Eq. (3.6), the moments of the spectral function, the r derivatives of the Euclidean correlator 
at (3/2, and the time derivatives of the real-time retarded correlator at t — are in one to 
one correspondence. While a short time expansion can never be used to rigorously extract 
transport coefficients, they have proved useful in in non-relativistic contexts [20, 21] 

For times which are short compared to the collision time it is reasonable to expect that 
the motion of heavy quarks is described by the free-streaming Boltzmann equation. Even 
in the interacting theory, the free streaming Boltzmann equation will describe the first time 
derivative of the retarded correlator or the first term in the Euclidean correlator, Eq. (3.6). 

Let us create an excess of heavy quarks, and subsequently study the diffusion of this excess 
at short times. This can be done by introducing a small chemical potential /i(x) = /io+<5/i(x) 
as in Section II. Then the thermal distribution function at an initial time t — is 



/o(x, p) 



U + f p (i ± /„; 



(3.7) 



with 1 , fp = l/(e^ p_w ^ T =Fl) . For short times the collision-less Boltzmann equation applies, 



h v 

dt 



/(x,p,t) = 0. 



The solution to this equation with the specified initial conditions is 

/(x,p,t) = / (x- V'P) • 
Then the fluctuation in the number density is 

d 3 p 



5N(x,t) 



(3.8) 



(3.9) 



(3.10) 



with 5/(x, p, t) = /(x, p, t) — f p . Then taking spatial Fourier transforms with k conjugate 
to x and substituting the distribution function, Eq. (3.9), we have 



SN(k, t) 



d 3 p 

(2tt) 3 



e- ik ^f p (l ± f p ) 



(3.11) 



Generally we will restrict ourselves to a heavy quark limit where there are well defined high and low 
frequency contributions. The discussion in this paragraph and the previous paragraph applies whenever 
the scale separation persists, and is therefore applicable to relativistic weakly coupled quarks . We will 
therefore generalize this paragraph to relativistic quarks with Bose-Einstein and Fermi-Dirac statistics. 
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For small times, we expand the exponential, and find 



5N(k,t) = 

with 



,2 



6fi(k), (3.12) 



^ = ^ = ?/(S^ 1± ^« (3 - 13) 
and 

T) = f^/(g« 1± «f- < 3 ' 14 > 

Thus, from Eq. (3.6), Eq. (2.7), and Eq. (3.12), we find 2 

G^(k,r)=T Xs (k) . (3.15) 

In the free theory, at k = there there are no corrections to this result and the Euclidean 
correlator is a constant. At finite k, the lattice correlator is not a constant even in the free 
theory. For massless particles, (v 2 /3) = 1/3, while for massive we have (v 2 /3) = T/M. 

We have outlined the short time expansion of XnnO*-, t) . Further insight is gained from 
the full free spectral function. From, Eq. (3.11), Eq. (2.7) and a simple Fourier transform 
we deduce that the retarded correlator from the free streaming Boltzmann equation is 

Xnn(K U) = l= f f p (l ± f p ) ^ . • 

T J (2ir) 6 uj — k ■ Vp + xe 

Taking the imaginary part, the corresponding spectral density is 

73 T 



_low 



1 f crp 

%u) = - J ^^/p(l±/p)k-Vp£(w-k-v p ). 



As shown in Appendix B, this form for the spectral density is identical to the one loop 
spectral function of the free theory at small k and lu, Eq. (BIO). As discussed in Appendix B, 
the resulting integral can be performed in the non-relativistic limit and we find the free 
spectral function for the heavy quark current-current correlator 

pfr(k,a;) = x.^ , 1 = exp[-— ^- I . (3.16) 



k 2 



This is the dynamic structure factor of a free non-relativistic gas [21]. In the free theory, 
the spectral function is essentially a Gaussian, with a width that is proportional to k 2 . In 
the limit that k = the correlator is 

p^(k,uj) = Xs( V 4-)^)- ( 3 - 17 ) 



2 Here we have considered only a single component gas. For the case of heavy quark diffusion, Eq. (3.13) 
and Eq. (3.14) should be multiplied by 4iV c to account for the sum over spin, color, and and anti-quarks. 
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In the free theory, the low frequency spectral density is infinitely narrow at k = 0. 
The moments of the spectral density are in one to one correspondence with the derivatives 
of the Euclidean correlator at (3/2. Since higher moments of a delta function are zero, 
all derivatives at (3/2 vanish and the low frequency contribution of the free theory to the 
Euclidean correlator is simply a flat line. Thus, provided the high frequency contribution of 
the spectral function can be subtracted, any bending of the Euclidean correlator is indicative 
of something beyond free streaming. In the next sections we will discuss how interactions 
smear the 5(u) function and estimate how much the Euclidean correlator curves at (3/2 as 
a function of diffusion coefficient. 



IV. HEAVY QUARK DIFFUSION IN THE LANGEVIN EFFECTIVE THEORY 

In this section we will discuss the predictions of the Langevin equations for the retarded 
correlator. As mentioned before, the time scale for heavy quark transport, M/T 2 is much 
larger than typical time scale for light degrees of freedom in the plasma. For this reason 
we will assume that the Langevin equations provide a good macroscopic description of the 
thermalization of charm quarks [12], 

dx l p 1 
~dt ~ ~M ' 

^ = at) - VP 1 , (e(t)e(t')) = Ks»6(t - o . 

The drag and fluctuation coefficients are related by the fluctuation dissipation relation 

For timescales which are much larger than 1 / 77 the heavy quark number density obeys 
ordinary diffusion equation 

d t N + DV 2 N = . 

The drag coefficient rj can be related to the diffusion coefficient through the Einstein relation 

Mr) k v ' 

The effective Langevin theory can be derived from kinetic theory in the weak coupling 
limit [12] and probably is adequate for describing heavy quark diffusion even for strongly 
interacting plasma. The Langevin equations make a definite prediction for the retarded 
correlator. Following the framework of linear response, consider an initial distribution of 
heavy quarks when a small perturbing chemical potential is applied, /i(x) = /io + <fyt(x). 
The initial phase space distribution of heavy quarks is 

/(x, p, t — 0) = e^'^e'Mff . (4.3) 
Summing over spins and colors, the initial number density of quarks minus anti-quarks is 

7V(x, t = 0) = [4N C ] — e-r sin h I J . (4.4) 
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By comparing Eq. (4.4) and Eq. (2.8), we find the static susceptibility 

X, = [4iV c ] (^) 3/2 ^cosh(^). (4.5) 

Let P(x, t) be the probability that a heavy quark starts at the origin at t — and moves 
a distance x over a time t . Consider the relaxation of an initial distribution of heavy quarks 
iV(x, t — 0) slightly perturbed from equilibrium. The distribution of heavy quarks at a later 
time is, 

A^/^Ptx-x',,) *(*<,«,), (4.6) 

or 

N(k,t) = P(k,i)JV(k,0). (4.7) 

Comparing this result with the linear response result, Eq. (2.10), we conclude that for small 
k and times large compared to typical medium timescale 

X™(k,t) = - Xa (k)$P(k,f). (4.8) 

Thus, to find the retarded correlator XnnQ*., uj) , we need only find the probability P(x, t). 

The probability distribution P(x, t) is determined in Appendix A. Not surprisingly, the 
distribution is a Gaussian, 

P(x ' t)= (Mt))3/ xp H^))' (4 ' 9) 

with a width that depends non-trivially on time 

2D 

<x 2 (t) =2Dt-— (l-e-*). (4.10) 

For large times, we have <r 2 (t) « 2Pt as expected from the ordinary diffusion equation. For 
small times, we have 

a 2 {t)^{T/M)t 2 (»7t (4.11) 

which reflects the initial thermal velocity distribution of heavy quarks, (v 2 /3) = T/M. Using 
Eq. (4.8), the probability distribution Eq. (4.9), and the definition of the retarded correlator, 
we find the following form: 

POD 

X JVJ v(k,cu) = Xs (k) / dte iwt k 2 D(l-e- r ' t )e- k2Dt + ^ D / r '^ 1 - e - rit K (4.12) 
Jo 

Eq. (4.12) summarizes the contribution of the Langevin equations to the retarded density- 
density correlator. The retarded correlator has following properties: 

1. For small k, Dk 2 <ti r], and arbitrarily large times, we may write the integrand as 
k 2 D (e~ k Dt — e~ nt ), and perform the integration 

Xs(k)Dk 2 Xa (k) Dk 2 

Xnn(K uj) = . — — — . 4.13 

— iu + k 2 D — %uj + i] 
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FIG. 1: The spectral density of the longitudinal current-current correlator 7rpjj(k, lu)/lu divided by 
.Dxs(k) as a function of a scaled frequency u) = uoD (M/T) for various values of a scaled momentum 
k = kD-y/ M/T. The solid lines show the spectral density from the Langevin equations for non- 
zero k. For comparison, the dotted lines show the spectral function of the free theory, Eq. (3.16), 
expressed in the same k and u> of the interacting theory. The dash-dotted line shows the k = 
result of the Langevin equations, Eq. (4.14). 



For small frequency u ~ Dk 2 , the first term dominates and recalls the diffusion equa- 
tion, (pt + Z?V 2 ) -1 . For large frequencies u> ~ r], xnn recalls the drag term of the 
Langevin equations, (dt + rf)" 1 . Of particular relevance to lattice measurements is the 
spectral density of the current-current correlator at k = 

pjj{0,u) _ 1 lmxjj{0,uj) _ T 1 rj 

UJ 7T UJ M 71 US 1 + r\ l 

2. The typical relaxation time of a heavy quark is set by the inverse drag coefficient, 
= D (M/T). The typical distance that a heavy quark moves over the relaxation is 
^/T/M/t] = Dy/M/T. The correlator 

Xnn is a function of a scaled spatial momentum 
k = \s.D\J M/T and a scaled frequency u = ooD (M/T). In Fig. 1 we show the spectral 
weight of the current-current correlator. For comparison we also show the free current- 
current correlator from Eq. (3.16). 
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3. Noting that x™(k,w) = ~ e iuJt x s d t P{k, t) with P(k,t) = e -k 2 ° 2 (t)/2^ it is easy to 
verify the consistency relation x(k, 0) = Xs(k)- 



V. NUMERICAL ESTIMATE OF THE EUCLIDEAN CORRELATOR 

In this section we will give a numerical estimate of the Euclidean vector current correlator. 
We will parametrize the spectral density with low and high frequency contributions. 

PjJ (k, co) = p l ™(k, co) + p*f(k, co) . (5.1) 

The high frequency part is present at zero temperature and will be parametrized as a 
J/ip resonance plus a continuum 



p-f (k = 0,.) = MlJl 5(co> - Ml,) + ^ - AMI) ^^-^ (| + A jjf) ■ 

(5.2) 

Here f v is the J/ip coupling to dileptons as described in Appendix C. The continuum 
contribution is motivated by the free spectral function calculated in Appendix B, but we 
have replaced 2M with the open charm threshold 2M D . 

For the low frequency part of the spectral function we will take two functional forms. 
The first form is the Lorentzian from the Langevin equations 

PjJ (k = 0,u>) T 1 V 

CO M 7T U0 Z + T) 1 

where r\ = This form is rigorously true when <C T, and the frequency small 

co<r]<^T. 

These inequalities are strained in our numerical work. For instance, for T/M c 1/5 and 
D ~ 0.25/T, is not really much less than T. Further, as discussed in Section III, the 
transport contribution is dominated by the second moment of the spectral function 



/ 



-PjA^co*. (5.4) 

to 



For the Lorentzian, this moment diverges and the transport contribution to the correlator 
is sensitive to the high frequency behavior of the ansatz where the Langevin approach is not 
valid. The higher moments open up the white noise in the Langevin equations. We therefore 
considered a Gaussian ansatz which falls much more rapidly at infinity 



(co) T ■ 



X s — e 2v g . (5.5) 

oo As M v /2^| ^ > 

The parameter, r] G = \f\^^i is fixed from the relation between the spectral density and 



the diffusion coefficient coefficient, 



_ xsD_ rj^g integral under this smeared delta 

function is again XsT/M. By comparing these functional forms we obtain a feeling for the 
uncertainties our the estimate. 
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FIG. 2: The ratio Gjj(t, T) /G j j (r, T) for different temperatures and k = 0. 



The temperature dependence of the Euclidean correlators comes from two sources: the 
temperature dependence of the spectral function p(k, cu,T), and the trivial temperature de- 
pendence of the integration kernel, Eq. (2.16). We obviously want to separate the interesting 
temperature dependence coming from the spectral function from the trivial temperature de- 
pendence coming from the integration kernel. This can be done by defining the reconstructed 
correlator [29]. 

r° cosher -1/(2T))) 

Gjj{Kr,T)= dupjj [k,w,T = . . 5.6 

Jo smh(w/(2T)) 

If the spectral function does not change above the deconfinement temperature T c , the ratio 
Gjj(k, r, T)/G r J e j(k, r, T) should be unity 

First we estimate the relative importance of the transport contribution to the correlator. 
For closer comparison with existing lattice data, we consider the diffusion of heavy quarks 
in a gluonic plasma where the transition temperature is T c = 270 MeV [32]. At this stage 
we can set the rj to zero (D = oo) and consider only the free spectral function. The charm 
quark mass M c is taken to be 1.3 GeV. In accord with lattice data [27-29], we will assume 
that J/ip is not modified by the medium and determine f v from its dilepton width (see 
Appendix C). M J/V , and Mr, are taken from the Particle Data Book [33] In Fig. 2 we show 
Gjj(k, t, T)/Gjj(k, r, T) for several temperatures. The transport contribution is of order 
7 — 12% and is the the only source of the temperature dependence seen in Fig. 2. A similar 
enhancement was found in actual lattice calculations [34]. 

Analytic understanding can be gained by performing the integral over the kernel at r = 
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1/(2T). In the heavy quark limit, we set M J/4 , w 2M and M D M, and find 

2 



^(k = 0,r,T) U/2 = 4iV c ^y /2 e ^^ + M*(f* 



2M 



M 



transport resonance 



, . MT\ ' m ( T 
4N r ~ e" ( 1 — — 



continuum 



f v /2M ~ 0.131 is small and suppresses the resonance contribution. The transport contri- 
bution is smaller by a factor of T/M relative to the continuum contribution. 

Interactions will modify the correlator by only a few percent. These small changes due to 
the transport must be disentangled from other in-medium effects such as a small shift in the 
mass or width of the resonance. This can be done by introducing a small chemical potential 
for the heavy quark, ji c <C M. Since the transport contribution is proportional to Xs, the 
small chemical potential will enhance the transport by factor of cosh(jU c /T), see Eq. (4.5). 
The small charm chemical potential will not affect the resonance and continuum contribu- 
tions to the spectral function to leading order in the heavy quark density, ~ e -( M -M/T 
Thus we expect that 

5Gjj = Gjj(t, T, /x) - Gjj(t, T, 0) (5.7) 

is largely insensitive to the high frequency behavior of the spectral function. For a thousand 
gauge configurations, the statistical error in the vector current correlators can be reduced 
below, 0.5%. One may hope that the same holds for the difference of the correlators, 5Gjj. 
Clearly, to achieve this precision one should difference the two correlators before averaging 
over gauge configurations. This needs to studied with numerical experiments. 

In Fig. 3(a) and (b) we show this difference for T = 1.1T C , /j, c = M/5 and different values 
of the diffusion constant D. As seen in Fig. 3, and as expected from Eq. (3.6), the effect 
of the diffusion coefficient is to provide a small curvature to the correlator and to shift the 
value of the correlator downward at r = (3/2. 

First we will concentrate on the curvature. If the final precision is 0.5% and D < l/(irT), 
then from Fig. 3(a), one could hope that the curvature is large enough to be determined 
in lattice simulations. In practice, it will be difficult to guarantee that the continuum 
contribution will not affect the extracted value. 

The downward shift of the correlator at (3/2 from the constant value, XsT/M, is a much 
larger effect. To isolate this transport contribution we consider the difference, 5Gjj(M), as 
a function of the heavy quark mass. We plot the ratio 



R{M) _ 5Gjj(M)/(xs(M)T/M) 



SGjj(M )/( Xs (M )T/M ) 



(5.9) 

=0/2 



For the free theory this quantity is one and is independent of the heavy quark mass. Devia- 
tions from one are a signature of interactions. Fig. 4(a) and (b) show this ratio as a function 
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FIG. 3: The difference of correlators at fi c = M c /5 and \i = for the (a) Lorentzian and (b) 
Gaussian ansatze and various values of the diffusion coefficient, D. 



of the heavy quark mass for the Lorentzian and Gaussian ansatze. Examining Fig. 4, we 
conclude that if the diffusion coefficient is sufficiently small, D < l/(nT), the transport peak 
should be visible in the mass dependence of the Euclidean correlator. Additional critical 
remarks are left to the conclusions. 



VI. BRIEF SUMMARY AND DISCUSSION 



The Euclidean current-current correlator is remarkably insensitive to the heavy quark dif- 
fusion coefficient. Indeed, to leading order in T/M, the Euclidean current-current correlator 
is independent of the diffusion coefficient. 3 This is explained as follows (see Section III). 



3 This is true whenever there is a separation between the transport and temperature time scales. Previously, 
Aarts and Resco found that Euclidean stress tensor correlations are independent of the coupling constant 
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FIG. 4: The relative transport contribution to the correlator as a function of quark mass, Eq. (5.9), 
for the (a) Lorentzian and (b) Gaussian ansatze to the low frequency spectral function. The 
numbers indicate the diffusion coefficient for each ansatz. 



The r derivatives of the euclidean current- current correlator at r = (3/2, the moments of 
the spectral function, and time derivatives the real-time retarded correlator at t = 0, are 
in one to one correspondence. Thus, the value of the current-current correlator (i.e. the 
zero-th derivative) is determined only by short times and may be calculated with the free 
streaming Boltzmann equation. In the end, the value of the current-current correlator at 
(3/2 is simply XsT/M, where Xs is the static susceptibility and T/M reflects average thermal 
velocity squared. Higher r derivatives (or moments of the spectral density) reflect the width 
of the transport peak and contain useful information about the transport time scales. 
In a free theory, the spectral density is proportional to a delta function 



PJj(k 



0,u) 



to 



which reflects the fact that in the free case, the spatial current is conserved in addition to 
the charge. This result may be found either by using the free streaming Boltzmann equation 
(see Section III) or performing a one loop expansion (see Appendix B). Since the spectral 
density is proportional to a delta function, higher r derivatives, or moments of the spectral 
function, vanish and the Euclidean current-current correlator is a constant, independent of 
r (see also Ref. [26]). In the interacting theory the delta function is smeared. Using the 
Langevin equations of motion, we analyze in Section IV how this delta function is smeared 
as a function of k and u. This result together with the free theory is summarized in Fig. 1. 
At k = 0, the Langevin effective theory dictates the replacement 

x( \ 1 ^ 



7T oj 2 + rf ' 



to leading order [26]. 



16 



where 77 = T/(MD) and D is the diffusion coefficient of the heavy quark. 

With this Lorentzian form for the spectral function at small omega, we adopted a simple 
transport + resonance + continuum model for the full spectral function and studied how 
the Euclidean correlator is modified by the transport peak in Section V. We also smeared 
the delta function with a Gaussian to illuminate the sensitivity to the Lorentzian ansatz 
which is only valid in a heavy quark limit and for u ^ T / (MD). 

Generally, the transport contribution to the full correlator is suppressed by a factor of 
T/M relative to the continuum contribution (see Eq. (5.7)). To disentangle the trans- 
port from the continuum and resonance contributions we proposed differencing two current- 
current correlators - one at finite heavy quark chemical potential and one at zero chemical 
potential, SGjj(t) = Gjj(r,fi) — Gjj(r,0) . This difference is proportional to the low fre- 
quency contribution and is independent of the high frequency contribution to leading order 
the heavy quark density, ~ e -( M -^)/ T _ With this procedure, the transport contribution can 
be separated from the other contributions at least parametrically. In practice (as opposed 
parametrics) our numerical work in Section V shows that extracting this piece is difficult 
though not impossible. A major unknown is the final precision when the difference of correla- 
tors is calculated. Clearly, one should difference and then average over gauge configurations. 
Exploratory lattice studies are needed to estimate this precision. 

The transport contribution to the correlator is displayed separately in Fig. 3(a) and (b) as 
a function of the diffusion coefficient. As analyzed in Section III, the effect of the diffusion 
coefficient is to shift the value of the current-current correlator down from its free value 
XsT/M, and to curve the correlator at (3/2. Parametrically, these effects are suppressed by 
(T/MD) 2 relative to XsT/M. The figure illustrates that if the diffusion coefficient is much 
greater than 1/T it will be difficult to measure the second derivative at (3/2. However if the 
precision is 0.5% it may be possible, although it will be hard to guarantee that the continuum 
contribution has been completely subtracted. To eliminate the continuum contribution it is 
desirable to make the mass as large as possible. On the other hand, the transport signal is 
proportional to (T/MD) 2 and therefore is suppressed by the mass. Ultimately, numerical 
experiments will determine the optimal heavy quark mass. 

Even with these complications, Fig. 3 shows that the Euclidean correlator at (3/2 is clearly 
shifted downward from its free value, XsT/M. This shift also is indicative of the width of 
the transport peak. To evaluate the magnitude of this shift, we proposed measuring 

SGjj(M)/( X s(M)T/M)\ T=m , 

as a function of quark mass; this quantity is independent of the mass in the free theory. 
As is shown in Fig. 4(a) and (b), in the interacting theory the width of the transport peak 
makes this quantity mass dependent. Judging from Fig. 4, if the diffusion coefficient is less 
than < 0.25 /T the effects of the transport peak should be visible in this mass dependence. 

Measuring Fig. 3 and Fig. 4 on the lattice is very difficult. The importance of such 
measurements should spur effort. Only measurements of this kind can seriously challenge 
the strong coupling assumptions that underly the hydrodynamic interpretation of the RHIC 
results. 
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APPENDIX A: DIFFUSION OF A BROWNIAN PARTICLE 



The goal of this appendix is to determine the probability P(x, t) that a Brownian particle 
will move a distance x from the origin over a time t. Consider the discretized Langevin 
equations: 

x t+1 - x t = |i , (Al) 
Pm - Pt = -VPt At + 6 , > = |^<W , ( A2) 

where the noise is drawn from a Gaussian distribution with the specified variance. 

Let W^[po, Pi, • • • , p n ] be the probability of having a sequence of momenta, po, Pi, • • • , Pn, 
where po is the momentum at time zero and p^ is the momentum after N time steps. The 
probability of having momentum p is given by the thermal distribution 



g 2MT 

P (Po) = ( 2vr MT) rf / 2 • 

Here and below d = 3 is the number of space dimensions. The probability to have momentum 
pi given p is the probability that the noise will attain the appropriate value 



At \ d/2 

P(pi|po) = J dttPfa- fa- ripote + 
Continuing in this way we deduce that probability distribution is 



W[p , Pl , • • , p n ] = {2nMT)d/2 {27rKM)Nd/2 exp \- + ™T) • ( A3 ) 

where p, = (p i+ i - pi) /At. 

Now the probability to move a distance Ax over a time At can be written as 

» N N-l 

P(Ax, At) = H d d Pi W[p , Pl , . . . , p N ] 5 d (Ax - J2 Af) . (A4) 

i=0 i=0 

We now rewrite the delta function as a Fourier integral and substitute Eq. (A3) into Eq. (A4) 
to obtain 

2 

JcZl n - P " 1 

d k TT^„ „ik.Ax e 2mt i ^ 



(2vrMT) d / 2 (2vr/tAt)^/ 2 



/ ,^At ^At , 2 \ 

x exp I -i 2^ — k • Pi - 2^(P« + VPi) • (A5) 

V i=0 i=0 / 

The integrals in Eq. (A5) are all Gaussian and can be performed. We performed the integrals 
in reverse order, p n , p n _i, . . . , pi, po and finally the k integral. The result is a Gaussian 
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with width 



r 2 _ T T 2 , K j 



where the discretized integrals I\ and I2 are, 

AT-l 



/1 = At ^(l-TjAf)' — > [ dt' 

AT-l 

/2 = (At) 3 £ 



-r,(t-t') 





to 


-> / dt' 

Jo 


fdte-^-^ 


-J=0 






Jf 



Performing the continuum integrals, and liberally using the relations D = = yields 
our final continuum form for the width: 



2D 

a 2 (t) =2Dt- ^(l-e-*) 



(A6) 



For large times, we have <J 2 (t) 2Dt as expected from the ordinary diffusion equation. For 
small times, we have cr 2 (t) ~ (T/M)t 2 reflecting the initial thermal distribution of heavy 
quarks, (v 2 /3) = T/M. 



APPENDIX B: THE FREE SPECTRAL FUNCTION 



To evaluate the high frequency behavior of the spectral function let us evaluate the free 
spectral function using standard methods [31]. To this end we will calculate Matsubara 
correlator 



, ft, 4 ; 



J dr j d 3 xe-^- k ' x (J£(x,r)J£(0,0)> , 



(Bl) 



with = k E = 2irnT. With this definition of the Matsubara propagator the real time 
retarded propagator can be determined from its Euclidean counter part through the relation 



X ^(k, A; ) = (-i) r G%(k, -ih -k° + i v ) , 



(B2) 



where r = <5 M o + 5 v q is the number of zeroes in the indices //, v. In the notation of the rest 
of the paper XAw(k, a;) = x°°(k, oj) and XjjQ*-, oj) = fc'fc- 7 'x' 3 '(k, cu). 

The one loop contribution to the spectral function is shown in Fig. 5. 



-j, + M) 

IE 



(#-t + M) 



p 2 + ei ' E (h 4 - P4 y + E 2 _ p 



Ye 



(B3) 



Here indices are raised and lowered with the metric tensor g E v = diag(-l,-l,-l,-l). 7^ satisfies 
{t£, Y E } = 2g% and j> = p^ E = -p^% - f lE - 
Let us examine a typical term in Eq. (B3) 



J„(k, -ik A )=TY,P2 



P4 



V\ + El (h ~ V,) 2 + E 2 p _ k ' 



(B4) 
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FIG. 5: Feynman graph contributing to the free spectral function 



where n — 0, 1, 2. Performing the frequency sum [31] we have, 

-l(+iE p ) n 



J n (k, -ifc 4 ) 



p^p-k 



1 - n p - n p _ k 
— ik,A — E n — 



-l) n (l - np -n 



p-kj 



+ 



p -"-^p-k 

(-i) n K 



-ifc 4 + -Ep + £p-k 

Wp-k) 



— i/c 4 - E u + E„ 



—ikd + E„- 



(B5) 



Evaluating the correlator in Eq. (B3) involves performing the trace, evaluating the frequency 
sums with Eq. (B5), and performing the continuation —ik^ — > k° + irj as indicated by 
Eq. (B2). The only contribution to the imaginary part of the correlator comes from energy 
denominators. In Eq. (B5) for example, the imaginary part of a typical energy denominator 
after the continuation — ik 4 — > k° + irj is 



Im- 



(k° + ir)) - E p - E p _ k 
With this identity we have 



ir5(k° -E p - £ p _ k ) . 



Imx 00 (k, k°) 



7T 



N r . 



lmx ij (k, k°) 



IT 



d 3 p 



(2vr) 3 4£ p £ p _ k 
d 3 p 



[(4E p k°)D + + (-8E 2 p + 4 p • k) D_] , (B6) 



Nr. 



(2tt)3 4£ p £ p _ k 



[(4^A; j + 4ibV - 8pV - 4p k 5 ij ) D 
+ {AE p k°)^DJ\ , 



(B7) 



where the even and odd functions D± are 

D ± = (1 - n p - n p _ k ) (S(k° -E p - £ p _ k ) ± 5(k° + E p + E p _ k ) 
+ (np - n p _ k ) - E p + £ p _ k ) ± + £ p - £ p _ k )) . 



(B8) 



The first pair delta functions can only be satisfied when \k°\ is large \k°\ ~ 2 M. The second 
pair of delta functions can be satisfied when \k°\ ~ k. Thus for k C T, the full correlator 
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can be written as a sum of high and low frequency contributions 



Imx^k, fc°) 



7T 



lmx^(k° k) 



7T 



+ 



low 



Imx<* u (k° k) 



7T 



high 



First let us focus on the high frequency contribution to the spectral density. To reach 
an analytic expression for the spectral density we set k = 0. Then the integral over 
(S(k° — 2E P ) ± 5(k° + 2E p )) is easily performed, yielding 



IT 



high 



llmx"(k = 0,o;) 
3 



7T 



high 



AW 
8tt 2 



2 ' 4M 2 /2 4M 2 , 

1 — I - + — — I tanh 



tar) ' 



(B9) 



This agrees with an earlier calculation [35] after accounting for a factor of two which results 
from a sum over two flavors in that calculation. 

Next we consider the low frequency contribution to the correlator which comes from 
difference of energies, 5(k° — E p + -E p _k)- For kcTwe expand to first order, 



dn 



dE r 



k - v p , 



with v p = p/Ep. Then the spectral density is 



Imx 00 (k, A; ) 



7T 



low 



d 3 p N c 



^3 ^ t V*° ("g) k • Vp - k ■ V P ) + * ( *° + k . Vp)] 

+ 8E 2 k • v P [S(k° - k ■ v p ) - 6(k° + k • v p )] } . 



Integrating over cos(0fc p ) eliminates the combination of delta functions symmetric with re- 
spect cos(#fcp). Integrating the anti-symmetric combination of delta functions yields a factor 
of two and therefore 



Imx 00 (k,u;) 



7T 



low 



_fp_ m ( dn 



dE r 



k • v p 5 (u - k • v p ) 



(BIO) 



Eq. (BIO) is identical with the correlator deduced from the free streaming Boltzmann equa- 
tion. 

This expression for the retarded correlator is readily simplified in the non-relativistic limit 
where n p = exp (— p 2 /(2MT)). The delta function can be written as 



k -v P %-k-v p ) = — 8[co s6 pk - —je[p-—j . 
Integrating Eq. (BIO) we find a Gaussian with a width that is proportional to k 2 , 



(Bll) 



Im X 00 (k,cu) 



7T 



: exp 



UJ 



low 



(B12) 
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Here, (f 2 /3) = T/M and Xs is the static susceptibility in the non-relativistic limit, Eq. (4.5). 
In the limit that k — > the width of the Gaussian approaches zero and we have 



Imx 00 (k = 0,uj) 



Ti 



(B13) 



low 



With this knowledge and the relation between the density-density and current-current cor- 
relators Eq. (2.22), we find Xjj 



Im xi/(k,w) 



TV 



Xs 



1 



low 



k 2 



2nk 2 



. cxp 



In the limit that k — > this function also approaches uS(u) 



Ti 



low 



(B14) 



(B15) 



APPENDIX C: RESONANCE SPECTRAL FUNCTION 



The coupling of a J/ip to the electromagnetic current at T = can be written as 

\p,a)=eQf v M J/ ,e»(p). (CI) 



Here M J/4 , is the J ftp mass, J^ M = eQ ay^c , e the charge of the positron, Q = +2/3 and f v 
is the electromagnetic decay constant. In writing this equation we have used the fact that 
Pn (0|^sm(0)Ip> a ) vanishes by current conservation. The decay rate of unpolarized J/tp into 
e + e~ may be expressed in terms of f v : 



r(J/V> -> e+e") = 



4tt Q 2 a 2 
3 M '""' /! 



EM r2 



(C2) 



Using the Particle Data Book [33] we obtain, f v /M J/4 , = 0.131. 

Using Eq. (2.12), Eq. (2.14), and Eq. (2.21), the spectral density at k = can be written 
as follows: 



where Df { (k, a;) is 



;(k,w) = J d 4 ie +iuMk ' x (j ! (x) J*(0)) 



(C3) 



(C4) 



Here the averages denote thermal averages and J^(x) = cy^c. We will assume that the J/^ 
coupling and mass are independent of temperature and simply replace the thermal average 
with vacuum averages. In the frequency domain of the resonance we may assume that one 
particle intermediate J/ip states dominate the correlator. Inserting one particle states we 
find 
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Using translation invariance, (0\J l (x)\pa) = e ip ' x (0| J*(0)|pcr), we perform the momentum 
and space-time integrals and find 

D ^^ = W k 5{u ~ Ek) E<°l Jt (°)l kff >< k ^l^(0)|0> • (C6) 
We now specialize to k = and use Eq. (CI) to obtain 

-6{u - M m ) f v M* . (C7) 



3 2M, J/i>iJv j/tp 



L J/4> 

A similar calculation yields D^(0, 00) and the resonance contribution to the spectral function 
reads 

f 2 M 2 , 

Pjj(0, u) = [S(u - M m ) - 5{u + M Jf ,)} . (C8) 
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